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In  this  final  report  we  summarize  our  conputational  efforts  addressing  the  thermal 
decomposition  of  3-nitro-l , 2 ,4-triazol-5-one  (NTO)  and  hexahydro-1 , 2 ,3-trinitro- 
1,3,5  triazine  (RDX) ,  These  efforts  expand  our  earlier  work  on  the  commone  explosive 
2,4,6-trinitrotoluene  (TNT),  which  has  been  described  in  Report  AL-TR-90-004  (32 
pages) .  We  also  report  on  the  dimethylnitramide  (DMNA)  decomposition,  bonding 
properties  of  bridged  hydrogens,  and  a  superacid  solvation  model.  Tautomer ism  is  j 
the  common  thread  in  the  studies  on  NTO  and  RDX,  which  are  summarized  in  sections  1 | 
and  2,  respectively.  In  section  6  we  expand  the  earlier  reported  imine  -  enamine,  \ 
keto  -  enol,  and  nitro  -  sci-nitro,  tautomerisms  with  the  nitroso  -  oxime  tautomerism 
using  GZ  theory  with  specific  inclusion  of  the  effects  of  solvent  molecules.  Isotope 
effects  are  often  used  as  indicator  of  the  decomposition  mechanism  for  the  expulsiort 
of  nitro  groups.  Hence,  we  report  the  douterium  kinetic  isotope  effects  for  decom- l 
position  of  dime thy lnitramine  using  denstry  functional  theory.  This  is  summarized  J 
in  section  3.  Superncid  initiated  electrophilic  nitrstions  reactions  are  intriguing 
processes,  which  are,  however,  not  modeled  conveniently  because  of  solvation  (pro- 


tonation) effects .  Therefore,  we  developed  a  computational  model  that  determines 
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the  superacid  solvation  effect  on  bases  and  electrophiles.  This  is  summarized 
in  section  5.  In  the  course  of  many  studies  we  became  interested  in  H-interaction 
energies  and  bonding  features.  We  decided  to  revisit  the  well  established  3 
center  ,  2  electron  interaction  system  which  underlies  many  non-classical  bonding 
arrangments.  This  fundamental  work  is  summarized  in  section  4.  Most  of  these 
studies  have  been  completed  successfully.  All  but  one  have  appeared  in  the 
literature.  These  will  tehrefore  discussed  only  in  a  summary  fashion.  A  more 
detailed  presentation  is  given  in  Appendix  A  for  the  nitroso  -  oxime  tautomerism. 
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i  INTRODUCTION 

i  In  ihis  Anal  report  we  summarize  our  computational  efforts  addressing  the  thermal  decomposition 
of  3-nirro-l  >?,4-triazol-5-one  (NTO)  and  hexahydro-l^^-trinitro-l^^  triazinc  (RDX),  These 
efforts  expand  our  earlier  work  on  the  common  explosive  2, 4, 6- trinitrotoluene  (TNT),  which  has 
been  described  in  Rcpun  AL-TR-90-004  (32  pages).  We  also  report  on  the  dimethyinitramide 
!  (DMNA)  decomposition,  bonding  properties  of  bridged  hydrogens,  and  a  superacid  solvation  model. 

[  Tautomcrism  is  die  common  thread  in  the  studies  on,  NTO  and  RDX,  which  are  summarized  in 
j  sections  1  and  2,  respectively.  In  section  6  we  expand  the  earlier  reported  inline  enanaine,  keto 
j  enol*  and  nitre  ••  aci-nitro,  tautomerisms  with  the  nitroso  *+  oxime  tautomcrism  using  02  theory 
|  with  specific  inclusion  of  die  effects  of  solvent  molecules. 

I  Isotope  effects  air.  often  used  as  indicator  of  the  decomposition  mechanism  for  the  expulsion  of  nitro 
|  groups.  Hence,  we  report  the  deuterium  kinetic  isotope  cITccls  fur  decomposition  of  dimethyl- 
!  nitraminc  using  density  functional  theory.  This  is  summarized  in  section  3. 

i  Superacid  initiated  electrophilic  nitrations  reactions  arc  intriguing  processes,  which  arc,  however, 
jnot  modeled  conveniently  because  of  solvation  (protonaiion)  effects.  Therefore,  we  developed  a 
i  computational  model  that  determines  the  superacid  solvation  effect  on  bases  and  electrophiles.  This 
*  is  summarized  in  section  S. 

i  In  the  course  of  many  studies  we  became  interested  in  H-interaction  energies  and  bonding  features. 
We  decided  to  revisit  the  well  established  3  center  2  electron  interaction  system  which  underlies 
|  many  non-classical  bonding  arrangements.  Tills  fundamental  work  is  summarized  In  section  4. 

i  Most  of  these  studies  have  been  completed  successfully.  All  but  one  have  appeared  in  the  literature* 

! These  will  therefore  discussed  only  in  a  summary  fashion.  A  more  detailed  presentation  is  given  in 
!  Appendix  A  for  the  nitroso  oxime  tautnmerisrn. 
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TECHNICAL  DISCUSSION 


1.  NT( 


Published  jlli  N.J.  Hams  and  K.  Lammertsma,  Tautomerism,  ionization,  and  bond  dissociations  of 
5  ~nitro-2/t-dihydro3H-l^,4-triazolone,  J.  Am.  Chem.  Soc.  118,  8048  -  8055  (1956). 


he  more  sal 


r>ninrs  are  snmr 


.  Tautomerization,  ionization,  and  bond  dissociations  of  the  insensitive  high-cncrgy  explosive  5*nixro* 
2,4-dihydro-3//-l,2,4-triazolonc  (NTO,  1)  were  studied  by  molecular  orbital  SQF  and  MP2  theories 
and  with  the  Recke  JLYP  hybrid  density  functional  using  the  6*3 1+G*  and  6-311+G**  basis  sets. 
Energies  computed  with  these  methods  were  compared  against  accurate  G2  energies  for  the 
tautomerization,  ionization,  and  bond  dissociations  of  nitromctliauc.  An  important  outcome  ftOtB 
|  the  present  study  is  our  finding  that  the  B3LYP  hybrid  density  functional  is  an  accurate,  and 
|  computationally  efficient  method  for  computing  energies  and  geometries  on  large  molecules.  Both 
MP2/6  31 1+G**  and  B3LYP/6-31 1 1 G**  gave  satisfactory  energies  for  laulouiciizatiun,  iunizitdon 
and  bond  dissociations  of  the  small  model  system,  nitromethane.  On  the  other  hand,  SO'  theory  is 
jnot  adequate  especially  for  computing  bond  dissociation  energies. 

jThe  B3LYP/6-31+<3+  and  MP2/6-31+G*  geometries  for  the  NTO  anion  9  compare  well  with  the 
j  reported  X-ray  structure  of  the  N 1X1  diaminoguanidinium  salt.  The  harmonic  viratioxuil  frequencies 
.calculated  with  Beckc3LYP  compare  well  with  those  observed  for  crystalline  NTO. 

1 

jThe  rather  small  energy  difference  between  NTO  and  the  hydroxy  tautomers  2  and  4  (4  keal/mol  at 
iMP2/6-311+G**  and  9  keal/mol  at  B3LYP/6-311+G**)  is  explained  to  result  from  enhanced 
| aromaticity  in  2  and  4  compared  with  1.  The  enhanced  aromaticity  of  these  l-H-1,2,4-  triazole 
!  systems  is  evident  from  calculated  C-N  bond  lengths,  chemical  shifts  of  the  NH  protons,  and 
imagnetical  susceptibility  anisotropies.  Likewise  the  rather  high  acidity  of  NTO  is  explained  in  terms 
i  of  the  delocalised,  at  ununic  ciuu  uucr  uf  die  NTO  anions.  The  computed  proton  affinity  for  the  N 1 X ) 
(anion  is  a  low  321  keal/mol  (B3LY P/6-3 1 1+G**  +  ZPE). 

(The  aci-nitro  tautomers  5  and  <3  arc  30  koal/mol  less  stable  than  NTO  itself.  The  difference  in 

(stabilities  is  considerably  larger  than  the  14  keal/mol  difference.  for  the  parent  tautomeric  system 

nitromethone/aci-nitromcthanc,  which  is  probably  a  consequence  of  repulsion  between  the  adjacent 

nitrogen  sp2  lone  pairs  in  5  and  6.  However,  die  uct-iiiuu  tautomers  are  Stable  enough  that  they 

jmight  play  a  pan  in  the  NTO  decomposition. 

1  . 

1  . 

jThe  estimated  N-II  and  C-N02  bond  dissociation  energies  fur  NTO  arc  respectively  93  and  70 
ikeal/mol,  and  that  for  the  N-OH  bond  of  the  tfei-nitro  tautomer  5  is  38  keal/mol.  Our  calculations 
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show  LliaL  several  tautomers  of  NTO  are  energetically  accessible,  and  these  may  play  a  role  In  the 
:ondensed  phase  process,  where  they  can  be  formed  by  a  base-catalyzed  mechanism.  Some  possible 
processes  in  the  initial  stages  of  NTO  decomposition  are  bimolccular  hydrogen  atom  transfers  (in 
he  condensed  phase),  C-N02  bond  hemolysis  (at  high  Leiupcmtmcs  or  under  conditions  of  shock 
yr  impact),  and  homolysis  of  the  N-OH  bond  in  the  nci-nitro  tautomers. 


Investigated  tautomers,  radicals,  en  anions  itfNTQ 


i 

! 
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2.  RDX 


Published  in:  N.J.  Hams  and  K.  Lanunertemo,  Ab  inilio  density  functional  computations  of 
conformations  and  bond  dissociation  energies  for  htMthydro-1 ,5  J-trlnltro- 1  JJ-mazene,  J.  Am. 
Chern.  SOC.  m.  65S3  -  6589  (1997). 

The  more  salient  points  arc  summarized; 

Ab  inido  computations  are  described  of  the  conformers  and  bond  dissociation  energies  of  hexahydxo- 
1.3,5-trinitro-1.3,5-triazlne  (RDX)  using  the  13  ecke3-Lee- Yang-Parr  (R3LYP)  nonlocal  density 
functional  and  the  standard  6-3 1G*  and  6-31 1G**  Gaussian  basis  sets.  The  B3LYP/6- 
31  lG**//B3LYP/6-3lG*  calculations  suggest  that  the  gas  phase  RDX  molecule  exists  as  a  mixture 
of  chair,  boat,  and  twist  conformers.  The  range  in  energy  for  these  conformers  is  only  about  t 
kcal/mol.  Th m  energy  minimum  chair  conformations  were  found,  differing  in  the  axial/equatorial 
orientations  of  the  nitro  groups.  Because  of  the  low  barrier  for  nitramine  inversions,  these  thrw 
conformers  should  interconvert  quite  rapidly ,  giving  on  average  chair  structure  with  three  equivalent 
(nitro  groups  and  C3v  symmetry.  This  is  in  agreement  with  the  experimentally  observed  gas  phase 

jstrucrure  for  RDX. 

\ 

[The  observed  conformation  for  crystalline  RDX  is  similar  to  the  conformation  for  the  theoretical 
[global  energy  minimum  structure.  Differences  between  observed  and  theoretical  N-N  bond  lengths 
jin  RDX  can  be  explained  by  crystal  packing  effects,  Intermoiecular  electrostatic  interactions  in  the. 
(solid  state  favor  resonance  donation  of  the  amine  lone  pair  into  the  NO^  group,  leading  to  shorter  N- 
|N  bonds  than  in  the  gas  phase.  The  RDX  molecule  hits  short  intramolecular  O—H  contacts,  which 
suggests  that  RDX  can  react  by  elimination  of  HONO. 

The  RDX  nitrogen  and  carbon  radicals  have  key  roles  hi  thermal  decomposition  of  RDX.  The 
B3LYP/6-311G**  computed  estimates  for  the  OH  and  the  N-NO2  dissociation  energies  (D0)  in 
RDX  are  respectively  85  and  42  kcal/mol.  The  latter  value  includes  a  correction  for  the  tendency 
Ifur  B3LYP/6-3 1 1G**  theory  to  underestimate  N-NO2  bond  strengths  by  6  kcal/mol  compared  with 
;liigh  level  G2  theory.  These  are  smaller  than  typical  secondary  OH  and  N-NO2  bond  strengths  in 
icyclohexane  (Vft  kcal/mol)  and  DMNA  (46  kcal/mol,  our  best  estimate).  Thus,  the  initiation  of 
'RDX  decomposition  by  N  NO2  homolysis  and  propagation  of  the  decomposition  by  hydrogen  alum 
abstractions  should  l>c  facile  processes, 

The  low  N-NO2  bond  strength  is  a  consequence  of  the  sigma  electron  withdrawing  NO2  substituents 
in  RDX,  The  weak  C-H  bond  in  RDX  is  similar  in  strength  to  an  allylic  or  benzylic  C-H  bond,  rattier 
than  a  typical  aliphatic  C-H  bond.  The  geometric  features  and  atomic  spin  densities  in  the  RDX 
jearbon  radical  radical  show  it  is  partly  delocalized  into  the  adjacent  NO2  groups. 
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Five  RDX  cortfarmers  with  R3LYP/6*$1G*  bond  lengths  (in  A)  and 
B3LYP/6-3J1G **  relative  energies  in  kcal/mvl. 
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3.  DMNA 

PuMishg^in;  N  J.  Harris  and  K.  ktmmertsma,  Ab  initio  density  functional  calculations  of  deuterium 
kinetic  isotope  effects  for  decomposition  ofdimcthylnitramine,  J.  Phys.  Client.  A  i01>  1 370  -  1373 
(1997). 


Dints  are  snmma 


Ab  inuio  harmonic  force  fields  were  computed  using  the  Becke3-Lee- Yang-Parr  hybrid  density 
functional  with  the  standard  6-3 1G*  basis  set  for  dimethylnimimine,  dimethyiamino  radical,  and  for 
the  transition  structure  for  five-center  elimination  of  HONO  from  dimcthylnitnuninc.  Hiesc  force 
fields  were  used  to  calculate  the  primary  deuterium  kinetic  isotope  effect  for  the  HONO  elimination, 
and  the  secondary  isotope  effect  for  the  N-N  bond  homolysis  ofdimcthylnitramine.  The  computed 
primary  effect  is  k’Vk136  =  4.21 ,  and  the  computed  secondary  effect  is  kn/k£>s  =  1 ,40,  both  at  240  *C. 
Comparison  with  the  experimentally  observed  isotope  effect  of  1.57  for  decomposition  of 
dimcthylnitramine-dO  in  solution  tu  240  "C  suggests  a  significant  part  of  the  observed  effect  i?  due 
to  the  secondary  effect  on  N-N  bond  homolysis.  Similar  beta  secondary  deuterium  isotope  effects 
are  expected  for  the  N-NOz  cleavages  that  initiate  the  decompositions  of  the  nili amine  explosives 
HMX  and  RDX.  However,  die  experimentally  observed  isotope  effects  on  the  global  decomposition 
|  rates  of  HMX  and  RDX  in  the  condensed  phase  are  best  interpreted  as  primary  effects  operating  on 
|  pathways  that  involve  rate-determining  C-H  (or  C  D)  homolysls. 


j  UdLYHIft-  Y  /  f?*  transition  structure  for  the  elimination  of  HONO  from  dimcthulnitramlnc . 

S 
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4. 3c-2e  Systems 


Published  in:  K.  Lammertsma  and  T.  Ohwarfn,  Three  center  -  two  electron  systems.  Origin  of  the 


lilting  of  its  substituent .v.,  J.  Am.  Chem.  Soc.  118. 7247  -  7254  (1996). 


jl  hc  more  salient  points  are  summarized* 


(This  computational  study,  using  MP2  6  31 1G**  theory,  is  concerned  with  species  that  display  the 
Acll-known  three  center  *  two  electron  bonding,  all  of  which  are  either  minima  or  transition 
Structures.  Particular  emphasis  is  pinned  on  the  origin  of  the  tilling  of  their  H-substituents.  Common 
bonding  patterns  are  found  for  a  broad  group  of  such  systems  that  include  mono*  and  dications, 
jnono-  and  dianiuns,  and  neutral  molecules  that  contain  eight  different  elements  of  the  first  and 
W.nnri  row  of  the.  periodic  table  The  parnnr  systems  studied  are.  XZH^  (C,  B,  Si,  Al),  (C, 

B,  Si,  Al),  and  XjH2*®2"  (C,  B,  Be.  Si,  Al,  Mg).  All  but  two  of  these  (i.e.,  the  ethyl  cation  and  the 
ethylene  diotuion)  have  their  hydrogens  bn  one  side  of  a  plane  that  contains  both  heavy  elements. 
Accordingly,  these  dements  have  inverted  geometries.  The  geometrical  tilt  angles  of  the  terminal 
lydrogenG  can  be  quite  significant  It  is  shown  that  the  geometrical  tilt  angles  between  the  X-II,  and 
lie  X-X  bonds  of  all  the  3e-2q  electron  systems  (anions,  cations,  dianions,  dications,  and  neutrals 
containing  eight  different  elements  of  the  first  and  second  row)  become  larger  when  the  3c-2e 
ystems  contain  the  heavier  elements.  Thus,  the  geometrical  tilt  angles  correlate  lineary  with  Allred 
Jfe  Rochow  s  elemental  electronegativities .  The  bond  lengths  between  the  heavy  elements  also 
correlate  with  the  tilt  angles.  The  origin  of  these  angles  is  traced  to  a  mixing  of  the  n  (X— X)  and  o 
XH(i))  orbitals.  These  two  orbitals  combined  show’  a  more  diffuse  nature  ill  the  H-bridguig  uicu. 
this  is  substantiated  by  the  one-electron  density  analysis,  which  shows  that  the  electronic  structures 
of  the  3c-2e  have  no  tilt  angles  between  the  X-Hj  and  X-X  bonds.  Instead,  molecular  graphs  show 
that  the  X-X  bonds  have  convex  curvatures.  They  further  show  that  the  bridging  hydrogens  have 
bond  paths  directly  to  the  centers  of  these  convex  X-X  bonds. 


A 


Group  A  Group  B 


-H 


Group  C 


CjHA.  B2Hs- 
SisHr>’.  AljHj' 


C2H3*.  B2Ha' 
SI2H,+.  AIjM,' 


C2H4~,  B2H4,  Be2H4" 
5l2H4+*tAl2H„MgzHft“ 


The  three  groups  of  this  study.  A:  XJIy~,  B:X2Hf'~,  and  C:  X2Hj 
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3,  Protonation 

Published  In:  NJ.  Harris,  T,  Ohwada.  and  K.  Lammertsma,  Hrainnaiion  energies  in  fluorosulfonic 
acid  using  ab  initio  seif  consistent  reaction  field  theory,  J.  Coinput.  Chcm.  ,12, 250  -  257  (1998). 

The  mure  salient  points  arc  summarized: 

Electrostatic  solvation  free,  energies  were  computed  for  several  small  neutral  bases  and  their 
conjugate  acids  using  a  continuum  solvation  model,  called  the  Self-Consistent  Isodcnsiiy  Polarizable 
Continuum  Model  (SCIPCM).  The  solvation  energies  were  computed  at  the  restricted  Hartree-Fnck 
(RIIF)  and  at  the  second-order  M0ller-P1esset  (MP2)  levels  of  theory,  as  well  as  with  die  Becke3- 
Lee-Yang-Parr  (B3LYP)  density  functional  theory,  using  the  standard  6-31G**  Gaussian  basis  set. 
The  RHF  solvation  energies  arc  similar  to  those  computed  at  the  correlated  MP2  and  B3LYP 
theoretical  levels. 

A  model  for  computing  protonotion  enthalpies  for  neutral  bases  in  fluorosulfuiuc  acid  solvent  leads 
to  the  equation  AHpro[tHS03F(B)  =  -PA(B)  +  AEl(BH1)-  AE,(B)  +  CONSTANT,  where  PA(R)  is  the 
ga.c  phase  pmtnn  nffmity  tor  base  B,  AEt(BH*)  is  the' SCIPCM  solvation  energy  for  the  conjugate 
acid,  and  AEt(B)  is  the  solvation  energy  for  the  base.  A  fit  to  experimental  values  of 
AHproi,HS03F<B)  fur  0:11  neutral  bases  (IIjO.  MeOH,  Me30,  H2S,  McSH,  MejS,  NH3,  MeNHU, 
Mc-NH,  and  PH9)  gives  constant  -  238-.4±2.9  keal/mol,  when  AAEj  is  computed  using  the  0.0004 
e_Bohr3  isodencity  surface  for  defining  the  solute  cavity,  at  the  RIIF/6-31Q++  level.  The 
comparison  with  available  experimental  data  for  aqueous  solvadon  Of  amines  and  ammonium  ions 
suggests  the  SCIPCM  energies  are  quite  reasonable,  especially  after  careering  experimental  da*?  for 
nonelectrostatic  contributions  to  the  solvation  energy.  Importantly,  the  SCIPCM  solvation  energies 
reproduce  the  expected  trend  of  decreased  solvation  encigy  fur  substitution  of  X-H  by  X-CH3.  OUT 
model  for  computing  enthalpies  for  protonation  of  neutral  bases  in  fluorosulfonic  acid  should  give 
enthalpies  accurate  to  within  5  keal/mol  for  second-row  bases  (thioe  there  and  phosphines)  and  2 
keal/mol  for  firet-row  bases  (amines,  ethers,  carbonyl  compounds). 

The  model  predicts  that  for  carbon  monbxide  AH^^^CO)  =  10  keal/mol.  Thus  protonation 
of  CO  is  endothermic,  nnd  the  conjugate  acid  HCO*  (formyl  cation)  behaves  as  a  Strong  acid  in 
fluorosulfonic  acid.  Thus,  the  model  allows  prediction  of  protonation  enthalpies  fa  other  bases  like 
carbon  monoxide  that  are  incompletely  protonated  in  superacids. 
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6.  Tautomerism 

Attachment  A:  J.A.  Long,  N.  J.  Harris,  and  K.  Lammeratna,  Formaldehyde  oxime  ni&osotnctharic 
tautomerism.  ;•  . 

The  more  salient  points  are  srnnma rized:  ' 

Formaldehyde  oxime  **  nitrosomcthanc  tautomerism,  isomeric  nitrone,  and  their  common  cations 
and  anions  were  studied  with  Gaussiaii-2  theory  using  MP2(full)/6-31G*  geometries  and  with 
density  functional  theory  using  B3LYP/6/31I+G**.  Geometrical  parameters,  harmonic  vibrational 
frequencies,  relative  stabilities,  conformation^  stabilities,  and  ionization  energies  are  compared  with 
experimental  gas  phase  data  when  available.  The  fuiiualdehydc  uxime  ^  nitrosomerhanc 
tautomerism  are  compared  with  the  amide  imidol,  imine  **  enaminc,  keto  **  end,  and  nitro  ^  aci- 
nitro  tautomeric  processes.  Solvent  effects  are  estimated  by  the  self  consistent  isodensity  polarizable 
continuum  model  (SCIPCM).  The  influence  of  hydrogen  bonding  interactions  with  the  solvent  is 
addressed  by  including  two  water  molecules.  In  the  final  evaluation,  fonnaldehyde  oxime  1  in  15.8 
kral/mol  more  stable  than  nitrosomethane  3  when  the  aqueous  solvation  correction  of  3.8  kcal/mQl 
ig  upplied  to  the  G2  energies.  Unsol vaied  farmoldchj'dc  oxime  is  estimated  to  be  11.1  kcal/mol  more 
stable  than  nitrone  4.  The  estimated  gas  phase  ionization  energies  (G2)  are  362.5  keal/mol  for 
formaldehyde  oxime,  350.6  kcal/mol  for  nitrosomethane,  and  351.4  kcal/mol  for  nitrone.  The 
corresponding  protonauon  energies  arc  191,6,  181,4,  and  202.8  kcal/mol,  respectively. 


Relative  energies  calculated  with  B3LYP/6-3 1G*  are  generally  in  very  good  agreement  with  those 
obtained  with  G2  theory.  Likewise,  their  geometries  and  harmonic  vibrational  frequencies  compare 


[well  with  those  obtained  with  MP2(fuliy6-3lG+. 
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